Theory of strong-field ionization of afigned CO2 



M. Abu-samha and L. B. Madsen 
Lundheck Foundation Theoretical Center for Quantum System Research, 
Department of Physics and Astronomy, Aarhus University, 8000 Aarhus C, Denmark. 

A theoretical framework for studying strong- field ionization of aligned molecules is presented, and 
alignment-dependent ionization yields are computed for CO2. Our calculations are in unprecedented 
agreement with recent experiments. We find that the ionization process is affected by intermediate 
resonance states, and the alignment-dependent ionization yields do not follow the electron density of 
the initial states. The theory explains the breakdown of semi-analytical theories, like the molecular 
tunneling theory and strong-field approximation, were excited electronic structure is neglected. 



Strong-field physics emerges as a very promising field 
for studying the structure and dynamics of molecular sys- 
tems. For example, recent developments in this field led 
to experiments on manipulating chemical reactions 
tomographic imaging of molecular orbitals [5], probing 
the nuclear dynamics on the attosecond and femtosec- 
ond [1] time scales, molecular alignment [5[ and torsional 
control 0. 

Molecules can be aligned relative to the laser field, and 
in order to fully exploit this effect, it is necessary to un- 
derstand the dependence of the initial ionization step on 
the molecular orientation. In 0, & B llS], alignment- 
dependent ionization yields were measured for ionization 
from the highest occupied molecular orbital (HOMO) 
of N2 K(2p)), O2 (^^(2p)), CO2 (^,(2p)) and CS2 
(7r^(3p)) in the transition regime between tunneling and 
multiphoton ionization, and a strong dependence on the 
alignment angle (/3 in Fig.[T]) was found. It was suggested 
that the angular dependence of the ionization directly 
maps to the orbital symmetry. The theoretical advance 
in the molecular case is impeded by the complexity of the 
molecular structure, which arises from the multicentre 
character of the molecular potential and the additional 
rovibrational degrees of freedom. Full ah initio calcu- 
lations of the alignment-dependent ionization are avail- 
able only for H+ pj, [jj, [jj, Uj] and H2 EE]. For 
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FIG. 1: (Color online) Illustration of the orientation (/3) of 
the two degenerate HOMO orbitals of CO2 with respect to the 
linear polarization axis {e) of the laser field. The contribution 
to the ionization yield from (b) is very small in comparison 
with that from (a), and is not considered in the present study. 
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FIG. 2: Radial molecular potentials ([T]) of CO2 up to /=8. 



larger molecules, despite a tremendous amount of exper- 
iments, no ah initio calculations are available, and the 
most widely used approaches to explain strong-field pro- 
cesses are the molecular tunneling theory [17] and strong- 
field approximation [l8|, ll9| . Calculations of alignment- 
dependent ionization yields based on these theories fail 
to explain recent experiments Q: Tunneling theory and 
strong- field approximation predict the ionization yield to 
follow the electron density of the initial electronic state, 
in contrast with observations for the CO2 molecule [9|. 

In this Letter, we use ah initio theory (within the 
single-active electron approximation) to investigate the 
response of linear, polyatomic molecules to intense fem- 
tosecond laser pulses, and in view of recent experimental 
results [9] we focus on the CO2 molecule. The time- 
dependent Schrodinger equation (TDSE) describing the 
active electron in the combined field of the frozen core 
and the laser pulse is solved numerically using grid meth- 
ods [20], and the time-dependent wavefunction is ana- 
lyzed using grid-based spectral methods [2l|, |22|] . With 
this theory, we are able to explain many features in the 
experiments [ol , and clearly identify the short-comings of 
the tunneling theory [l^ and the strong-field approxima- 
tion [isl,!!^. In short, the present work underlies the im- 
portance of dynamics in the excited state manifold, and 
extends studies within the single-active electron model to 
systems beyond H2", H2 and their isotopes. 
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FIG. 3: Ionization yields as a function of the angle p (see 
Fig. [T]) for CO2. The dashed (solid) line denotes calcu- 
lations based on the present approach for a 10-cycle laser 
pulse at 800 nm and laser peak intensities 5.6x lO^^W/cm^ 
(1.1 X lO^^W/cm^). The dotted line denotes calculations 
based on the molecular tunneling theory 0]. The computed 
yields are given relative to the orientation that gives the max- 
imum yield. 

The time-dependent wavefunction is solved in a partial 
wave expansion. Thus, it will be convenient if the poten- 
tial describing the active electron {V{r)) is also expanded 
in partial waves. V(f) is obtained from quantum mechan- 
ical calculations [23], and electron exchange is treated 
within the local-density approximation 0. The radial 
components of the potential, Uim{r)^ are obtained by in- 
tegrating over the angular variables, i.e. 

uur) = I dnvr^mvi^- (1) 

The integral is evaluated numerically using cubature on 
the sphere [25]. The radial potentials are shown for CO2 
in Fig. [21 The potential provides accurate description 
of the HOMO orbital (both orbital energy and angular 
decomposition) and the excited states of the molecule; 
the energies of the first and second excited states with 
TTu symmetry (5.85 eV and 2.72 eV) are in agreement 
with the values (4.75 eV and 2.48 eV) reported in x-ray 
absorption spectroscopy [26] (see also Fig. [5|). 

The TDSE is solved on a grid within the veloc- 
ity gauge |20]. The external field, linearly polarized 
along is characterised by the vector potential A{t) = 
Aq sin^ {irt / T) cos{Ljt)e with T the pulse duration, co the 
frequency and Aq = Eq/co with Eq the field strength. We 
use an equidistant grid with 2048 points that extends up 
to 160 au. The laser pulses contain 10 cycles, and the 
calculations were performed at 800 nm and peak inten- 
sity 5.6xlO^^W/cm^ and l.lx lO^^W/cm^. The angular 
basis set contains 21 spherical harmonics. The ioniza- 
tion yields are calculated from the numerical grid cal- 
culations by applying an absorbing boundary [21]. The 
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FIG. 4: Initial- state electron density (thick, solid line; axis 
to the left) and ionization yields (triangles; axis to the 
right) in strong- field ionization of CO 2 at laser intensities 
l.lxlO^'^W/cm^ (open triangles) and at 5.6x lO^^W/cm^ 
(filled triangles). 
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FIG. 5: Bound-state energy spectrum of CO 2 as computed 
at the peak intensity of a 10-cycle, 800-nm pulse (A(t)=0) 
with intensity (a) 5.6x lO^^W/cm^ and (b) l.lx lO^^W/cm^ 
The solid, dotted and dashed lines, respectively, denotes [3 
values 0°, 45° and 90°. In (a), the spectral structure from 
left to right is assigned to the following electronic states: tt^ 
(HOMO-1), Tig (HOMO), Tiu (LUMO) and ttu (C 3p). 

calculations were repeated in a larger box (320 au with 
4096 grid points) with a larger angular basis (31 spherical 
harmonics) and the results are converged. 

The ionization yields are shown in Fig. [3] as ob- 
tained from our TDSE calculations at 5.6x lO^^W/cm^ 
and l.lxlO^^W/cm^, and from molecular tunneling the- 
ory [9]. Starting with the TDSE results, the yields are 
very sharp and peak at 45±3°. The effect of laser inten- 
sity is rather negligible; it only enhances the ionization 
yield somewhat when the molecule is aligned parallel to 
the laser field. The molecular tunneling theory predicts 
the ionization yield to peak at /5=25°. 

In Fig. m we compare the ionization yields at orien- 



tation to the electron density of the HOMO orbital 
at that orientation. One clearly sees that the ionization 
yields do not follow the electron density of the initial 
state, e.g., the HOMO orbital density for CO2 peaks at 
/3=20° whereas the ionization yield peaks at 45±3°. This 
indicates that excited electronic states, accounted for in 
the present approach, contribute to the ionization pro- 
cess. To investigate this point further, we computed the 
bound states energy spectrum 2l|, [22[ 
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where w(t) is a hanning function and C{t) the auto- 
correlation function of the time-dependent wavefunction 
The bound-state energy spectra are shown for 
CO2 in Fig. [3 as obtained from the calculations at 
5.6xlO^^W/cm2 and l.lx lO^^W/cm^. The energy spec- 
tra are similar for the two laser intensities; the spectra 
show that ionization occurs via several excited states, 
namely the lowest unoccupied molecular orbital (LUMO) 
at energy 0.21 au below the threshold, and the carbon 
7r^i(3p) excited state at energy 0.10 au below the thresh- 
old. The excited states were identified by comparison 
to theoretical oxygen K-edge x-ray absorption spectra of 
CO2 [2^. The population in the LUMO state depends 
strongly on /3, and for P values up to 80°, it is significantly 
larger than the ionization yields. At 5.6x lO^^W/cm^, 
the total ionization yields (LUMO population at the 
peak of the pulse; A(t)=0) at /3 values of 0°, 45°, 
80° and 90°, respectively, are 7.35x10"^ (2.01 xlO-^), 
2.58x10-^ (1.34x10-2), 2.43x10-"^ (1.33x10"^) and 
5.31x10"^ (9.78x10"'^). The energy spectrum shows 
that some of the HOMO population gets transferred into 
the HOMO-1, 7r^^(2p) orbital. A crude estimate of the 
relative contribution from this state to the total ioniza- 
tion yield is about 10%, based on a simple model that 
takes into account the relative population and the differ- 
ence between the ionization potentials of the HOMO and 
HOMO-1 orbitals. 

Now we compare our calculations with the experimen- 
tal measurements for CO2 at l.lx lO^^W/cm^ [9]. First, 
it should be emphasized that while the present calcula- 
tions are carried out at a fixed laser intensity, the mea- 
surements take into account focal volume effects, fiuctu- 
ations of laser intensity and alignment distribution. In- 
cluding focal volume effects and fiuctuation of laser in- 
tensity places heavy demands on the computational re- 
sources, and since the ionization yields show little depen- 
dence on laser intensity (cf. Fig. [3]), these contributions 
will be left out. However, we do convolute our results 
with the alignment distribution function devised in [9|. 
The final results are shown in Fig. [6l We also apply the 
same convolution procedure to the ionization yields ob- 
tained from the tunneling theory [9|]. One can clearly 
see that the tunneling theory is no where near the mea- 
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FIG. 6: (Color online) Experimental (circles) and theoretical 
(dashed line: 5.6x lO^^W/cm^; solid line: 1.1 x lO^^W/cm^; 
crosses: MO-ADK theory) [9] ionization yields for CO2. 



surements. By contrast, our TDSE calculations show ex- 
cellent agreement with the measurements. This demon- 
strates the important role of excited states in strong-field 
ionization of CO2. 

In the preceding paragraphs, we demonstrated superi- 
ority of the proposed approach over the molecular tun- 
neling theory and the strong-field approximation. This 
is because the present approach takes into account the 
molecular electronic potential, and hence is capable of 
describing the ionization process in different ionization 
regimes. The tunneling theory [17] and strong-field ap- 
proximation [18, 19], on the other hand, assume a pure 
ionization mechanism (pure tunneling or direct multipho- 
ton ionization), without accounting for excited electronic 
structure. Therefore, they fail to explain the experimen- 
tal measurements. 

While the discussion above focused on the alignment- 
dependent ionization only, the present methodology can 
be extended to extract other observables of interest. For 
example, the above-threshold-ionization spectrum can be 
determined by extending the analysis based on the auto- 
correlation function to the continuum [2l| , the photoelec- 
tron distribution can be determined by fiux analysis \2?\ 
and the full momentum distribution can be determined 
by projection on exact scattering states jll|, or by anal- 
ysis based on asymptotic projection operators [28]. Also, 
since the expectation of the dipole acceleration is read- 
ily evaluated from the time-dependent wave packet, the 
method can be extended to consider the process of high- 
harmonic generation in molecules. 

In conclusion, our theoretical calculations are in un- 
precedented agreement with recent experiments, and ex- 
plain the breakdown of the tunneling theory and strong- 
field approximation. The main finding is that the excited 
electronic structure is vital in strong-field ionization of 
molecules, and the measured ionization yields do not fol- 
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low the electron density of the initial state. The dynam- 
ics in the excited state manifold will affect not only the 
ionization dynamics but surely also the process of high- 
harmonic generation and in particular the mapping to the 
electronic continuum of the excited state symmetry may 
impede the extension of tomography techniques to more 
general classes of molecules. The proposed framework for 
treating polyatomic molecules is grid based, which is the 
most widely used approach in strong- field physics, and 
takes input potentials from standard quantum chemistry 
codes [23|. It could therefore be easily implemented by 
many research groups, and will facilitate studies of sys- 
tems more complex than H2 and H2" and their isotopes. 
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